CO 

o 

C\! 

Cl 

o 

•<3 

Gk 


NRL  Report  7930 


A Program  to  Plot  an  Annotated  Track 
or  a Track  and  Bathymetry  or  Magnetic  Profile 
on  a Mercator  Projection 

Marilyn  L.  Blodgett  and  James  V.  Massingill 

Propagation  Branch 
Acoustics  Division 


February  27,  1976 


NAVAL  RESEARCH  LABORATORY 
Washingiost,  D.C. 


Approved  for  public  release;  distribution  unlimited. 


security  classification  of  This  pace  fl«i»n  Lmim  Enimd) 


h REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3 RECIPIENT'S  CATALOG  NUMBER 

Tj 

“A  TRACK  AND  BATHYMEriTY  OR  MAG^/ItIC 
^OFILE  ON  A ilERCATOR  PROJECTION  « 


6.  PERFORMING  CIR^  REPORT  NUMBER 


Marilyn  L.^lodgett  James  Vy^assingill  I 


9 PERFORMING  ORGANIZATION  NAME  AND  AOORESS 


Naval  Research  Laboratory 
Washington,  D.C.  20375 


1 1 I CONTROLLING  OFFICE  NAME  ANO  AOORESS 
Naval  Material  Command 
Via  Arctic  Submarine  Command 
Naval  Undersea  Center 


10  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a WORK  UNIT  NUMBERS 

NRL  Problem  SOl-47 
Project  ZF  52-552-001;  62759N 
^pa  Order  1787 


Feba 


>76 


& 


San  Dieeo.  CA  92132 

!'•  IWNI  wSmCTfb^CY  NAME  a AODR?SSfrni7BJ>lffiwifW»»I»W*BW«BWV»)/ 


’ CLASS,  (oi  thl»  report) 

Unfefafesified 


15*.  DECLASSIFICATION/OOWNCRAOING 
SCHEDULE 


If  nuiPiiptmwp^  TTyriiLTilWHIlimiJ  ' } 


Approved  for  public  release;  distribution  unlimited. 


t7  OlSTRieuTlON  STATEMENT  (oi  tho  obotfct  onUrod  in  Block  20,  U dUl^rotU  from  RopoH) 


18  SUPPLEMENTARY  NOTES 


It  KEY  WORDS  fCon(/nu«  on  rovotto  tido  If  nocoBOory  md  tdontliy  by  block  numbor) 


Oceanographic  data 
Geophysics 
Profiles 
Bathymetry 


Magnetic  surveys 
Navigation 
Mercator  projection 


ACT  ^ConUnuo  on  rovoroo  §ldB  It  nocBBBsry  mid  tdontlfy  by  block  mmbor) 

W program  has  been  written  for  plotting  an  annotated  track  or  for  plotting  a track  and 
the  superimposed  bathymetry  or  magnetic  profile  on  a Mercator  projection.  The  program 
reads  the  data  (navigation,  bathymetry,  or  magnetics)  from  a magnetic  tape  in  BCD  form. 

'^The  program  will  annotate  every  point  or  every  nth  point.  Navigation  is  annotated  with 
fix  numbers,  bathymetry  will  uncorrected  fathoms,  meters,  or  corrected  meters,  and  magnetics 
with  the  residual  magnetic  intensity.  _ ^ 


JConUnue^ 


DD  ,:Ty>  1473 


EDITION  OF  I NOV  6>  IS  OBSOLETE 
S/N  0 102-014-  6601  I 


1 SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whm  Dtim  BntwAdJ 


-/  96~C  ^ 


CONTENTS 


1.0  IDENTIFICATION  1 

2.0  PURPOSE  2 

3.0  USAGE  3 

4.0  METHOD  OF  ALGORITHM  5 

5.0  SOURCE  LANGUAGE  LISTING  5 

6.0  COMPARISON  5 

7.0  TEST  METHOD  AND  RESULTS  5 

8.0  REMARKS  5 

APPENDIX  A — Sample  Input  Data  Records  6 

APPENDIX  B - Deck  Assembly  9 

APPENDIX  C- Sample  Profiles  15 

APPENDIX  D - Flow  Chart  16 

APPENDIX  E — Source  Language  Listing  17 


iii 


PRECEDING  PikGE  BLANK-IWT  FILMED 


A PROGRAM  TO  PLOT  AN  ANNOTATED  TRACK 
OR  A TRACK  AND  BATHYMETRY 
OR  MAGNETIC  PROFILE  ON  A 
MERCATOR  PROJECTION 


1.0  IDENTIFICATION 

1.1  Title 

A Program  to  Plot  an  Annotated  Track  or  a Track  and  Bathymetry  or  Magnetic 
Profile  on  a Mercator  Projection. 

1.2  Identification  Name 
Mercator. 

1.3  Classification  Code 
None. 

1 . 4 RCC  Identification  Num_b_CT 
None. 

1.5  Entry  Points 
MERCATOR. 

1.6  ftogramming  Language 

Language:  CDC  3600/3800  Fortran. 

Routine  Type:  Program. 

Operating  System:  Drum  Scope  2.1. 

1.7  Compu^  and  Configuration 
CDC  3800. 

1.8  Contributor  or  Programmer 

Marilyn  L.  Blodgett,  Code  8176MB,  Propagation  Branch,  written  for  the 
Environmental  Sciences  Section,  Acoustics  Division. 

1.9  Contributing  Organization 

NRL  — Naval  Research  Laboratory,  Washington,  D.C.  20375. 

1.10  Program  Availability 

If  supplied  with  a magnetic  tape,  the  Environmental  Sciences  Section, 

Acoustics  Division,  will  make  a copy  of  this  program. 


Manuscript  submitted  September  3,  1975. 
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1.11  Verification 

This  program  has  been  used  and  tested  by  the  Environmental  Sciences  Section, 
Acoustics  Division,  for  several  months. 

1.12  Date 

May  1976. 


2.0  PURPOSE 

2.1  Description  of  the  Routine 

This  program  reads  the  data  collected  by  an  oceanographic  or  geophysical 
experiment  from  a magnetic  tape  and  plots  an  annotated  track  or  plots  the 
track  and  bathymetric  or  magnetic  value  perpendicular  to  the  track  as  a 
profile.  We  use  tlie  format  recommended  by  the  National  Research  Council 
of  the  National  Academy  of  Sciences  with  one  slight  modification  for  the 
input  data  tape.  There  is  one  logical  record  (of  80  characters)  for  each  data 
point.  The  different  types  of  data  (navigation,  bathymetry,  and  magnetics) 
are  separated  by  an  end-of-file  mark  with  a double  end-of-file  mark  at  the 
end  of  all  the  data. 

Before  the  program  reads  this  input  tape,  it  reads  three  cards.  The  first  card 
specifies  the  minimum  and  maximum  latitude  and  longitude  in  degrees  and 
minutes.  The  second  card  specifies  the  kind  of  plot  desired,  the  number  of 
tick  marks  between  the  meridians  and  parallels,  the  spacing  Detween  the 
meridians  and  the  parallels,  the  height  of  the  map  to  be  drawn,  the  dates  of 
the  data  on  the  first  input  tape  to  be  considered  for  plotting,  the  actual  values 
to  be  plotted,  and  the  units  per  inch  for  plotting  the  bathymetric  or  magnetic 
profiles  along  the  track.  The  third  card  defines  the  actual  data  format  on 
the  input  tape  (the  format  varies  for  the  three  types  of  data). 

With  all  the  required  parameters  defined,  the  program  starts  to  read  the  input 
tape  one  record  at  a time.  Each  record  is  checked  to  see  that  the  fix  falls  on 
the  defined  grid  and  that  it  was  taken  on  or  between  the  two  specified  dates. 
Only  those  points  which  meet  both  requirements  are  plotted.  The  program 
continues  reading  the  input  tape  until  it  reads  an  end-of-file  mark  or  a fix  taken 
after  the  last  specified  date.  If  there  are  additional  input  tapes,  the  program 
reads  them  in  a similar  manner.  The  beginning  and  end  dates  for  each  new 
input  tape  are  required  on  an  Extra  card.  A maximum  of  four  input  tapes  can 
be  used.  When  all  the  input  tapes  have  been  read,  the  program  prepares  to 
plot  the  track  or  annotated  track  and,  if  required,  prepares  to  plot  the  profiling 
values  for  either  bathymetry  or  magnetics. 

The  track  is  plotted  on  a Mercator  projection  which  is  drawn  exactly  to 
scale.  The  grid  may  be  blown  up  to  any  reasonable  size.  Since  the  projection 
is  drawn  exactly  to  scale,  a mosaic  can  later  be  built  of  the  entire  area. 


2 


NRL  REPORT  7930 


2.1.1  Navigation  Data 

The  program  reads  the  date  (month  and  day),  time  (hour  and  whole  minutes), 
latitude,  longitude,  and  fix  number  from  the  input  tape  according  to  the 
specified  fonnat.  The  southern  latitudes  and  the  western  longitude  are  preceded 
by  a negative  sign.  Normally,  the  track  is  plotted  in  a continuous  line  with 
every  nth  fix  marked  with  a square  symbol  and  annotated  with  the  fix  number. 

2.1.2  Bathymetry  Data 

The  program  reads  the  date  (month  and  day),  time  (hour  and  whole  minutes), 
latitude,  longitude,  and  uncorrected  fathoms  or  corrected  meters  from  the  input 
tape  according  to  the  specified  format.  The  southern  latitudes  and  the  western 
longitudes  are  preceded  by  a negative  sign.  The  program  can  convert  uncorrected 
fathoms  to  uncorrected  meters.  If  an  annotated  track  is  desired,  each  fix  is  plot- 
ted with  a small  plus  symbol  and  annotated  with  uncorrected  fathoms,  meters,  or 
corrected  meters.  Tf  a profile  is  desired,  the  track  is  plotted  in  a continuous 
straight  line,  and  the  profiling  series  is  multiplied  by  -1  to  drop  it  below  the  track. 

2.1.3  Magnetic  Data 

The  program  reads  the  date  (month  and  day),  time  (hour,  and  whole  minutes), 
latitude,  longitude,  and  residual  magnetic  intensity  from  the  input  tape  accord- 
ing to  the  specified  fonnat.  The  southern  latitudes  and  the  western  longitudes 
are  preceded  by  a negative  sign.  If  an  annotated  track  is  desired,  each  fix  is 
plotted  with  a small  plus  symbol  and  annotated  with  the  residual  magnetic 
intensity.  If  a profile  is  desired,  the  track  is  plotted  in  a continuous  straight 
line,  and  the  profiling  series  is  residual  magnetic  intensity. 

2.2  Problem  Background 

This  program  allows  the  researcher  to  plot  vast  numbers  of  bathymetric  and/or 
magnetic  data  on  a Mercator  chart.  There  are  two  ways  of  plotting  the  data: 
plot  an  annotated  track  (data  value  plotted  according  to  geographical  coordi- 
nates) or  plot  a profile  of  the  data  using  the  track  as  a reference  point. 


3.0  USAGE 

3.1  Calling  Sequence  or  Operation  Procedure 
Not  applicable. 

3.2  Arguments,  Parameters,  and/or  Initial  Conditions 
Not  applicable. 

3.3  Space  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage: 

7074  octeil  (3644  decimal)  locations  exclusive  of  system  library 
functions. 

3.3.2  Common  Blocks: 

/MAPSYS/,  /INPT/,  /RATCOM/,  /CHTLBL/. 
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3.4 

3.5 

3.6 

3.7 

3.8 

3.9 

3.10 

3.11 

3.12 

3.13 

3.14 


3.3.3  Temporary  Storage: 

None. 

Messages  and  Instructions  to  the  Operator 
None. 

Error  Return,  Messages,  and  Codes 
Insufficient  information  to  complete  plot. 

Informative  Messages  to  the  User 
None. 

Input 

The  actual  format  of  the  data  on  the  input  tape,  the  kind  of  plot  desired,  and 
the  map  specifications  are  read  in  via  input  cards.  The  track  and  the  actual 
data  to  be  annotated  or  profiled  is  read  in  via  magnetic  tape  on  logical  units 
10  through  13.  Appendix  A presents  samples  of  our  data  formats  on  the  in- 
put tape.  Appendix  B is  a complete  description  of  the  input  deck  setup. 

Output 

The  program  prints  on  the  standard  printer  (logical  unit  61)  the  data  format, 
chart  parameters,  and  the  actual  number  of  points  plotted  for  the  track  and 
profile.  It  writes  the  plotting  instructions  on  a magnetic  tape  (logical  unit 
40).  Appendix  C presents  sample  profiles. 

Formats 

Appendix  B describes  the  program  deck  structure. 

External  Routines  and  Symbols 

SKIPFILE,  PLOT,  SQRTF,  SINF,  LOGF,  COSF,  TANF,  ATANF,  ACOSF, 
NUMBER,  SYMBOL,  XMODF. 

Timing 

The  time  required  depends  on  the  size  of  the  grid  and  the  number  of  data 
read  and  plotted. 

Accuracy 

The  grid  is  reproduced  exactly  to  scale. 

Cautions  to  Users 
None. 

Program  Deck  Structure 

Appendix  B describes  the  program  deck  structure. 
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3.15  References  — Literature 

N.  Maddage  and  S.J.  Kozloski,  Program  Scachart,  Naval  Research  Laboratory, 
Washington,  D.C.  M.L.  Blodgett  < nd  J.V.  Massingill,  “A  Program  for  Storing 
Oceanographic  Data  on  Magnetic  Tape,”  NRL  Report  7861,  March  1975. 

4.0  METHOD  OF  ALGORITHM 
Not  applicable. 

5.0  FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 
The  flow  chart  and  listing  are  given  in  Appendixes  D and  E. 

6.0  COMPARISON 

No  other  known  programs  are  available  for  comparison. 

7.0  TEST  METHOD  AND  RESULTS 

The  program  has  been  used  and  tested  successfully  on  a Calcomp  Plotter. 

8.0  REMARKS 
None. 
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APPENDIX  B 

Deck  Assembly  for  Program  Mercator 


'The  program  uses  a scratch  tape  on  logical  unit  05,  but  no  Equip  card  is  required,  since  the  drum 
is  used, 

the  Fortran  source  deck  is  used  instead  of  the  binary  deck,  a Fortran  card  is  required  after  the 
Equip  card.  !n  addition,  a Scope  card  and  Load  card  must  follow  the  source  deck. 
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Number 

Card  Title 

Column 

Number 

Description 

1 

Job 

1-21 

7/9  JOB,  Charge  No.,  ID  No.,  time. 
See  page  2-2  of  the  3600/3800 
Computer  System  Drum  Scope 
Manual. 

2-6 

Equip 

1-18 

7/9  EQUIP,  40=**,  WO,  LO 
7/9  EQUIP,  10=**,  RO,  HI 
7/9  EQUIP,  11=**,  RO,  HI 
7/9  EQUIP,  12=**,  RO,  HI 
7/9  EQUIP,  13=**,  RO,  HI 
10,  11,  12,  13,  40  = logical  unit 
numbers. 

RO  = read  only.  LO  = low  density. 
WO  = write  only.  HI  = high  density. 
See  page  2-3  of  the  3600/3800 
Computer  System  Drum  Scope 
Manual. 

7 

Program 

Deck  of 
Mercator 

This  is  the  main  program  with  asso- 
ciated subroutines.  If  the  Fortran 
source  deck  is  used  instead  of  the 
binary  deck,  a Fortran  card  is  re- 
quired after  the  Equip  card.  The 
Fortran  card  reads  7/9  FTN,  L,  R, 
X.  In  addition  a Scope  card  with 
SCOPE  starting  in  Column  10  and 
a Load  card  must  follow  the  source 
deck. 

8 

Run 

1-13 

7/9  RUN,  T,  P,  R,  M,  D 
T = time  limit  in  minutes. 

P = maximum  number  of  print  or 

write  operations  ‘ 

R,M,D  may  be  left  blank.  See  page 
2-15  of  the  3600/3800  Computer 
System  Drum  Scope  Manual. 

9 

Chart 

1-5 

28.0 

Degree  portion  of  southernmost 
latitude  for  the  chart  (28°30'N). 
(If  the  latitude  were  28°30'S,  a 
minus  sign  would  precede  28.0.) 

6-10 

30.0 

Minute  portion  of  southernmost 
latitude  (28°30'N). 
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Column 

Number  Card  Title  Number 

11-15 

16-20 

21-26 

27-30 

31-36 

37-40 

10  Parameter  1-5 

7-8 

9-10 


Description 

35.0 

Degree  portion  of  northernmost 
latitude.  (A  minus  sign  would  in- 
dicate the  latitude  is  south.) 

0.0 

Minute  portion  of  northernmost 
latitude.  (Both  degrees  and  minutes 
must  be  specified  in  decimal  form 
even  if  a zero  value  is  assigned.) 

-45.0 

Degree  portion  of  westernmost 
longitude  for  the  chart  (45°10'W).. 
(The  minus  sign  indicates  the  lon- 
gitude is  west.) 

10.0 

Minute  portion  of  westernmost 
longitude. 

-40.0 

Degree  portion  of  easternmost  lon- 
gitude (40°25'W).. 

25.0 

Minute  portion  of  easternmost 
longitude. 

4000  (right  justified). 

Units  per  inch  for  plotting  values 
along  the  track.  In  this  example 
a value  of  4000  gammas  would  be 
plotted  1 inch  above  the  track.  The 
remainder  of  the  profiling  data 
would  be  scaled  accordingly. 

60 

Number  of  tick  marks  between  the 
meridians. 

40 

Number  of  tick  marks  between  the 
parallels. 
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Number 


Column 

Card  Title  Number 

12 


14 


16 


19 


20-24 


26-30 


31-35 


Description 

0,  1,  or  2 

0 = every  fix  will  be  plotted  with 
a + symbol  and  annotated. 

1 = the  track  will  be  plotted  in  a 
continuous  line  and  every  nth  fix 
(columns  20-24)  will  be  marked 
with  a small  square  and  annotated. 

2 = plot  value  series  as  a profile 
perpendicular  to  the  track. 

0 or  1 

0 = use  value  as  read  from  the  input 
tape. 

1 = read  uncorrected  fathoms  from 
the  input  tape  and  convert  to  un- 
corrected meters. 

0 or  1 

0 = do  not  multiply  the  profiling 
series  by  -1. 

1 = multiply  the  profiling  series  by 
-1  to  drop  it  below  the  track. 

0,  1,  or  2 

Number  of  files  to  be  skipped  over 
on  the  first  input  tape.  There  are 
a maximum  of  three  files  on  our 
Geodata  tapes. 

12  (right  justified). 

This  value  designates  the  n of  the 
nth  fix  to  be  annotated.  This  value 
is  specified  only  if  there  is  a 1 in 
column  12. 


5.0 

Degree  portion  of  the  interval  be- 
tween the  meridians. 

0.0 

Minute  portion  of  the  interval  be- 
tween the  meridians.  (Both  degrees 
and  minutes  must  be  specified  in 
decimal  form  even  if  a zero  value  is 
assigned.)  In  this  example  the  spac- 
ing will  be  at  5-degree  intervals. 
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Number 


11 


Column 

Card  Title  Number 

36-40 


41-45 


56-60 


Description 

10.0 

Degree  portion  of  the  interval  be- 
tween the  parallels. 

0.0 

Minute  portion  of  the  interval  be- 
tween the  parallels.  In  this  example 
the  spacing  will  be  at  10-degree 
intervals. 

9.00 

Physical  height  of  the  chart  to  be 
drawn. 


64  1,  2,  3,  or  4 

Number  of  input  tapes,  with  the 
maximum  being  four. 

65-72  02250330 

Date  and  time  of  the  first  data  point 
to  be  read  and  plotted  from  the  first 
input  tape.  Columns  65  and  66  are 
for  the  month,  columns  67  and  68 
are  for  the  day,  columns  69  and  70 
are  for  the  hour,  and  columns  71  and 
72  are  for  the  minutes. 


73-80  02280830 

Date  and  time  of  the  last  data  point 
to  be  read  and  plotted  from  the  first 
input  tape.  All  data  between  the 
date  and  time  of  the  first  data  point 
and  of  the  last  data  point  will  be 
read  and  plotted  if  it  falls  on  the 
defined  chart. 

Format  1-?  (15X14,  1X14,  1XF8.4,  F9.4, 

10XF5.1)  This  format  should  be 
replaced  by  the  desired  input  format. 
The  format  must  be  enclosed  in  pa- 
rentheses and  left-justified.  Via  this 
format  the  program  reads  the  date 
(month  and  day),  time  (hour  and 
whole  minutes),  latitude,  longitude, 
and  value  to  be  annotated  or  pro- 
filed (fix  number  of  navigation,  un- 
corrected fathoms  or  corrected  meters 
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Number  Card  Title 


12  Extra 


Column 

Number  Description 

for  bathymetry,  and  residual  mag- 
netic intensity  for  magnetics).  The 
formats  for  reading  the  three  data 
types  on  our  input  tapes  are: 

• Navigation 

(15X14,  1X14,  1XF8.4, 
F9.4,  9XF5) 

• Bathymetry  - uncorrected 
fathoms 

(15X14,  1X14,  1XF8.4, 
F9.4,  10XF5.1) 

• Magnetics 

(15X14,  1X14,  1XF8.4, 
F9.4,  28XF5). 

1-4  0,  1,  or  2 

Number  of  files  to  be  skipped  over 
on  the  second  input  tape.  There 
must  be  an  Extra  card  for  each  ad- 
ditional input  tape.  Since  there  is 
a maximum  of  four  input  tapes, 
the  maximum  number  of  Extra 
cards  is  three. 


5-12 


13-20 


13  End  of  File 


02250330 

Date  and  time  of  the  first  data  point 
to  be  read  and  plotted  from  the 
second  input  tape.  The  dates  for 
the  first  input  tape  are  on  the  Map 
Parameter  Card. 

02280830 

Date  and  time  of  the  last  data  point 
to  be  read  and  plotted  from  the 
second  input  tape. 

Terminates  the  run. 
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APPENDIX  E 
Source  Language  Listing 


C 


C 


•c 


c 


c 


IQI 

to 


PR0GRAH  heritor 
REAL  LAT 

real  LATHINflATHAX'.LeNQHlN.LQNQMAX 
CIMENS10N  IBLr<25A)|AN0R<2)»HEAD(2),S(2^ 
DIMENSION  IPMUoJ 
C0MM0N/MAPSVJ/ 


MPCi 


CeRRH)« 

RR» 

del. 

ACHT, 

IPR| 

XMIN, 

YMIN, 

xhax. 

YMAXi 

MP# 

ICR, 

lA. 

RAH|» 

RAHI 

r.  alatnmi 

ALATNS* 

ALATXD,  ALATXM, 

ALATXS, 

I'l  AL0GNM| 

ALGGNSi 

AL6GXD,  AL0CXH, 

ALOQXS, 

K0MPi 

WIDTH, 

HITE,  DEGDVD, 

DEGDVH, 

1'.  NTIC, 

TlCLGf, 

ALTMD(5),ALTMM{5) 

,ALTMS<9) 

ALNMC(5),ALNMM{5),ALNMS<5)  alTSD<5),ALTSM(5),ALTSS(5j, 
AlNSC(5)»AUnSM(5)iALNSS{5),arI5),  VR(5) 


C0MM0N/flATC9K/NAVSYSi 

N0RATE, 

RATIDI5), 

CDLY(5), 

RATINCI3), 

1 

ARATI?),  PHIMC5), 

ALAMI5), 

PHlS(5), 

ALASI5), 

KNTHAX, 

2 

KNT,  BpS, 

NR, 

DELT, 

0SLN(5) 

C0MMON/CHTLBI./MINMER, 

MNMSR(90)»MAXMER» 

HINPAR, 

MNPARISO), 

1 

naxpar. 

NC, 

NNCJ50), 

nt:, 

NLY, 

2 

NY<5o!» 

NRY, 

XWl3o>* 

YK(5q), 

KXAXIS- 

9 

KYAXIS. 

lYN, 

xtib, 

YTI8; 

4 

UN 

C0MM0N/flXQRAT/PCL,N?IX 


IP  • 0 

inRST»0 

KSUNTfl 


A , 

PR  1 «l 


800 

9qo 

1000 

1100 

1200 

1300 

1400 

ifoo 

l6oo 

1^00 

1«00 


2100 

2200 

2300 


2400 

2500 


2$00 

27oo 

2®00 


63200 

63300 


ICR  * «0  83400 

iEXiO 

imo 

IUPg3 

ALATNSaO.O 

ALATXSaO.O 

ALeQN6»0,0 

AL0QXS(<O,O 

CEOODS*0,0 

CEGOV6fO,0 

Pl»3,l4t59 

K0MP»t 

NAVSYS«t 

NSRATErQ 

lYNfO 

LYN«l 

MPCit 

READ  IN  RANGE  6P  CHART  63600 

REAO(6O|1O1)4LATNO.ALATNH,ALATXD|ALATXH,AL6OND.AL0ONH|ALOGXDi 
tALeOXM 

r0RKA.(4r5,i*r6,i,r4,i»r6,iiP4,i) 

REAP! iCR, 001 )OlNCH,NT|C,NTtX',ILINE,|USE,lHULT,{TAPR|tPriDEODVOiOeGDvH,DEQD 
?CVM, 
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1 


1 


.lCEQDDD,DEGDDt',HrTE,N0TPiJUDYl.lTHl,jUDY2.lTM2 
601  r0RHAT(r5ilXi3lS.I3,ll>ilXF5,li3F5,l,lOXF5,2,!5l4) 

Ce  777  JX?iilTApE 
777  CALL  SK1PFILE<IT) 
ce  TO(20i300JIA 
20  RRboiV966329?66 

^S  B a 

C«B*«  ceNvERT  degrees,  minutes  and  SRC0NDS  TQ  RADIANS 
call  0T0R(ALATKDiALATNH,ALATNS,C3RRm) 

CALL  DT0R{ALATXDiALATXM,ALATXS,CORR(2J) 

CALL  DTOR(AL'EQND,ALeGNM,AL0GNS,CflRR(3l» 

CALL  OTORIAUeGXD,AL0GX«,AL0GXS,C0ftR<An 
call  DTGRinBCDVD,DEGDvH,nE00VS,0eu 
CALL  D7QR(06EODD,OFQDPH,DEQOOS,OOL) 
call  PLeTS(lEuF,25A,4g,29) 
call  PLOTIo,(;,g,fli3) 

Cm**  KERCAT  IS  THc  SuER9UT1NE  NAME  FQR  MERCAT0R  PROJRCTI0N 
too  C*l-L  HERCAT 

tF(WlDTH.LE,C,C,«R. HITE, UE, 0,0)00  T3  900 
700  CALL  GRATie 

IF«NTIC)719,713,710 
7ig  CALL  TICX 
715  |F<NTIX)26,3<,717 

717  CALL  TICY 
26  CALL  CHLA8EL 
25  MBl 

AMlNBOtO 

SECsO,0 

minusraooooocooooooooob 


XLAM«A8SF{ALATNM) 
K6Y»ALATN0,AAD, MINUS 
lF<K0Y,EO,MpAu5r27,28 

27  |,ATHtN»ALATNC*XLAM/60iO 
ce  T0  29 

28  LA’iMIN»ALATNC  ♦ XLAH/60,0 

29  XLAM»ABSF(ALATXM) 
KeY«ALATXD,AAO, MINUS 
ir(K0Y.EO,MliAUS)3O,3t 

30  LATmAX«ALATXCqXLAH/6o,0 
09  T9  32 

31  LATmAXbALATXC  ♦ XLAM/60,0 

32  XLAM*ABSF|AL'CGNM) 
K0Y»AL0GND,A''O, MINUS 
IF(K0Y,EO,HIAUS)  33,34 


33  L0NGM|N*AL9Q^D'•XLAM/6^,o 
00  T0  35 

L0NGM|NaAL0GAD  9 XLAH/60,0 
XLAM»ABSFIAL'CGXM) 
K0VsAL0r,XD,AAD, MINUS 
IF(X0Y,EQ,MP''US)36',37 
LONGMAX8AL0G>O*XLAM/6oiO 
00  T0  38 

LeNQMAXsAL9Q>D*  XLAM/60,0 
REAO(60,1000n7M 
1000  F0RhAT<ioA9) 

WRlTE(6i,555MrM 


34 

39 


36 

37 

38 


599 

tOlO 

toll 


FCRMATI*  DATA  PSRhaT  ,,,*i10A8) 
V.RITE(6t,10lC) 

FORMAT (IHCil^HCHART  PARAMETERS) 
HRITE(6i,i01t)LATHINiLATMAX 
FORMATIiH  ,2cHS3LTHMaST  LATITUDE 

I*  2 ^ 


•F1O.2.1OX2OHN0RTHM08T  LATITUDE 


3200 


75300 

75400 


75500 

73600 

75700 

759flo 


37oo 

3000 

3900 

6200 


I 
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J0i2  r0RHAT«lH  ,2CRWBSTH0ST  LONOITUPE  .f 1q ,2.1oX2oHEASTH85T  LQNGJTUDH 
1 «^10|2) 

KH«9 

I0P«2 

7A00  REAO(IT,IKH)lDAY»lTlMt#RlAT|RU0NQiAN0HAL 
IF(l0CHECK,n)7<OO,7O5l 
7c&l  lF(E0f »lT)ioC2il001 
7Aq1  lUPij3 

G0  T0  7A00 

1001  lFUDAY,LT,Jl.oyi)  G0  T0  7401 
|rUDAY,GT,d(,DY2)  G0  T")  10C2 
IF(IDAY,F.(5,0I  08  Ta  1002 

IFUCAY.EC.JLDYl.ANO.lYlME.LT.ITHl)  G0  t0  7401 
IF{;DAV,E!.,J(.DY2,AN0,jTjME,GT,ITM2)Qa  TQ  1002 
lF(RLAT,LT.tiTMIN)  G0  T0  7401 
IF<RLAT,GT,L*TMA)I)  go  T0  7A01 
IF(RL?NG,!.T»L0NQHIN}  00  T0  7401 
IF(RL0NG,GT A0NGPAX)  00  T0  7401 
ITRAKfilTRAK  ♦ 1 

CALL  P-0R(RL*T  lAMlN.SEC.RRLAT) 

CALL  OT0R(RL'CNG  ,AMlNt?EC,RRL0NG) 

CALL  M6RCTR(CiAOKT,RR,RRLATiRRL0NG,XDIM,YDIM) 

IF(IU3E,EG,i»6fl9«6je 
6(j9  ANGHAL^ANeMAL*t|82B8 
610  IFUMULT,EO,n«Xli6o8 
6*1  AN0mAl*,AN0MAL 
60e  IF(JFJRsT)30C1, 8000,8001 
8g00  XLASTsXplM 
YLASTiYDlH 
${1)?XDJH 
S{2>»YDIN 
tFiRSTn 

CALL  PL0T(XOJM»YClH,3j 
8001  *N0M(iJbAN0M<2> 

AN0m(2)*ANOHAL 
:F(IFIRST,E0,1J  503,5Q4 
5o3  XDIFFbXD!K-S(i) 

VDIFF|YD!H.S{2) 

1FIRST»2 
00  T0  19 

904  XOIFFbCXOIFF  * |X0Im»S<1>»)/2 
YDIFF*(YniFF  ♦ nOlM  * S(2)))/2 
19  READ<2)?AYAN2 (YDiFr.XDiFF) 

IF<HEAD(2)  ,L1 ,0>t‘6A0<2)g2»Pl*HeAO(2) 

|F(o,,LE,HEAE{2|)4u,408 
411  |F(HEAD{2|,LE|Pr/2;i07i40e 

4fl7  ISIONSI 
00  T0  23 

408  |F(3i>P|/2,,LT.HEA0(2n412i410 


412 

409 

410 
23 


6q9 

7oO 


|r(HEAD|2),Lt,3|«Pl)<09|4l0 

ISIQNbI 

00  T0  23 

ISIGN|!>1 

X0NE3S<l)*lSlGN«ANaM<l)*SIN{HEAD(2) J/OINCH 
XTHG«S(2)*ISIGN*ANOM(1)»C9S<HEAO(2))/OINCH 
1F<ILINE,E0,’2>  605^606 

KRlTE(o5,7oo)  X0NE,xTH0|lUP 
F0RMAT(2Fior<iU> 
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8<1)»XDJM 

S(2)»YDJM 

»<EADU>gHEAD(2) 

CAUL  ANGANSODIMiXUASTiYOIMiYUASTiANG) 

CALL  nXANfl(XOtMiYDlH,ANG|AN0HAU,KOUNT,lLjNE,lPT,IUP) 
K0UNT5K0UNT  ♦ 1 
XL*ST||X01M 
YUASTgYDIN 
60  T0  7400 
1002  lEXslEX  « 1 
REWIND  IT 

IP<N0TP,EO,JEX)  iC0  TO  999 

IT»IT  * 1 

READ<ICR,997)ITAPE,jUpYi,tTMi#JU0Y2,lTM2 

997  rORKAT(5I4) 

CO  998  ixil.lTApe 

998  CAUL  SKlprn-ElIT) 

Kh«5 

I0Ps2 

IUP*3 

60  T0  7400 


999  IP<lLlME,gO>S)407,800 

8q7  rewind 

1003  READ(g5,7OO)X6NEiXTW0|IUP 
ir(IOCHECX,oJ)ioo3',5fl9 
5o8  ir(Eor,o5)aoc*aoi 

8ol  1AN0M5IAN0M  ♦ J 

IP(|PIRST,E0,2»  305-,906 
?09  call  PLOT(X0AEiXTW0i3| 
|PIRST*3 


906 

800 


1020 


CALL  P1.0T(X0KeiXTH0,IUP) 

60  TO  Ioq3 
CALL  PL0TS<OiO> 

CALL  STrpPLOT 
«RITC{61,102C)ITRAK 
rORMAT(lH0il5HPR«QRAH  PLOTTEDi 
IFilLlNE.NE.S)  oe  TO  I0l5 
WRlT6(6l,lflJl)lAN0M 


l9,5Xl9HrtXES  roR  THE  TRACK) 


1021  F0RMAT<iHoi13HPRCORAM  PL0TTEuh8,9X22HP0INTS  TOR 
1015  5T0P 


THE  PROriLE) 


END 
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length 

ENTRY  POINTS  MERCATOR 

el;ck  nahus 

MARSYS 

INRT 

RATCSH 

CHTlOL 

rpcRAT 

EXTERNAL  3YMQ0L5 


Q8CENTRY 

THENO, 

OliCOAjOO 

08CST0PS 

08CDICT, 

SKjPPtLE 

DTOR 

PieT5 

puer 

HERCAT 

ORATIg 

TfCx 

TfCY 

CKLAQgL 

HBP.CTR 

ANCANO 


OOASg  SYM00LS 


PPAN0 

8T0PPl;«r 

ATAN2 

SPNF 

C0$r 

08C{FESP 

O0CfPl«C 


RE^, 

T3R, 

srt-; 

SL'«, 

5U 


GNSJNQL 


01743 

00869 


IDENT 


00021 

00133 

00064 

OO473 

00002 


MERCATOR 
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S'jB«OijTJNE  OMTie 


• T>IS  SU-i'^OUTINE  STiRcs  u*LUeS  !N  ARRAY  TPAT  ARE  USED  F0R  P1.0TTINQ 

* SiRCr.R  GRATJC'.JLE  TGR  A K?:RCA-*e«  AND  1‘jvERSE  mERCAT'^R  PRQJECTIOM 


« Ct-AHT, 


I' jCuHENTATIG'J  ?0R  EXPLANATIQN  0F  INVERSE  HERCATQR, 


cgmmon/mapsvs/ 


JLSGI 

t-pc, 


* >L^MC(5),AL.VHM(5),AUNMS<5j,Ai.r3D(: 

5 ^U'^SC<S),AI.NSM(5),ALNSS<9),aR(5>, 

CSM.NS.N/'FlXQRn/DCUiNTlX 


VR<5» 


100 

110 


120 


C 

C : : a : : ::  X : = : = 3 s : ■ s 3 s : 1 s s a a s s s s a a a BQRQER  e r 3 s 3 B * X : = = = a : a 3 s s X ■ X a e X ; = 1 3 
C 

CNEDEG  ? 3,m9926539/180,0 
CALL  PL0T(XMlN)Y;'lN,3j 
CALL  PLoT(XhAX|Y|-1»u2) 

CALL  PLgT,XMAx,Yf'AX-,2j 
CALL  PL0t!xpIn,YRAX,2| 

CALL  PL0T(XMlfi,YPlN,2J 
tr  <Di:L)9oo»9coij,(Jo 
G0  T0(lioiil(|9lo)MP 
L«  s 5 
KL  a 1 
NTRY  R al 

lF(l)EL,LT,aNEDGG)NTRy  R 1 
TEMP  a C0RRILK)/CEL 
ITEllP  a TEMP  + OiOOOl 
CIFF  A ABSPnEMP-lTEHp) 

AINCR  a DEL 

IF<Dlf'F,LT,DELA2C,O)C0  TO  i'JO 
IF(C0RR(LK)iLT,O,O)Gtl  T0  130 
AINCR  a 0GL*a,0»DIFFj 
CO  TO  140 
A|NCR  a D6L«tlFF 
ALCHK  a CSRR(L’<*l)^AlNCR/4,0 
GO  TO(150il«C)KL 
ALONG  a C?RR (3 ) lA INCR 
ISB  a 1 
GO  TO  170 

ALAT  8 C0RRU)*A:NCR 
I SB  8 3 
NP  8 3 
ISW  P X 
KSW  ■ 1 


130 

l4n 

190 


160 

r/0 


60000 

6Q100 

60200 

60300 

60400 

60500 

60600 

6o7oo 

60800 

60900 


C0RRUJi 

RQ- 

del* 

ACHT, 

61000 

IPR| 

XHIN, 

YMIN, 

XMAX, 

61100 

YH*Xi 
RAHi  1 

MP» 

RAHI 

ICRi 

61200 

6l300 

61400 

, ALATNMi 

ALATNSt 

ALATXDi  AI.ATXM, 

ALATXS, 

61S(/0 

■,  AL£GNH, 

ALOGNSi 

ALOGXOi  ALCGXP, 

ALOGXS, 

61600 

KOMP, 

NIOTH, 

HITE,  DEcrn'D, 

DEGOVM, 

61700 

. NTJC, 

TICLGTi 

A1.7MD(5),ALTMH(3) 

,AI.TMS(5), 

81000 

6l900 

62000 

62300 

62400 


63o00 

63200 

63300 

63400 

63500 

63700 

63800 

63900 

64000 

64100 

64200 

1<300 

64400 

64500 

646Q0 

64700 

64800 

64900 

63000 

63100 

63200 

63300 

63400 
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1«0 


t90 


GO  T0<iaOii9C)Ki 

CAM,  H1:«CTR(c*^CI<T,RR,C6PR(ISB),A10NO|XX.YY) 
CALL  PL0T<XXiYY|NP> 
jS«  * Np.i 
CO  T0  (183,2C0)JS« 
l83  CALL  gHFGLR{H0NOiXXiNTRY) 
along  » AL0!NG*OEL 
IF<AL9NQ.LT,alCHK)G0  T8  200 
LK  1 i 
KL  ? X 
MRY  f 2 
GO  TO  220 

CALL  MGRCTR{C,A0hT‘,RR,ALAT»C8RRnS8J»XX-,YY> 
CALL  PL0T(XX»YY|NP) 

JSW  » Np«i 
G0  T0(195,2OC)JSk 
195  CALL  CHF0Ln(>LATiYY|NTRY> 

AL*T  8 ALAT<<I:6L 
1F(ALAT,GT.ALCHK)Q0  T0  900 
200  1F(NP,NE,3)QIs  T0  210 
tS8  q ISB*K9W 
KSH  • 'IqKSW 
210  ISH  • 

NP  t NP*tSW 
00  T0<lflOil9O)Kl, 

STDELsOEL 
CELqOQL 

|F(DEL)900,9CO>t20 
CALL  PL0T(XHlNiYPlN,3j 
CSLiSTDEL 
return 

2N0 


22q 


900 

«10 


65500 


658fln 

65900 

66000 

66100 

66200 

66300 

66400 

66500 


66900 

67000 

67100 

67200 

67300 

67400 

67500 

67600 

67700 

67800 

67900 


71400 


PfiSGRAP  LENQTH 

OO32I 

ENTRY  P0JNT5 
8LGCX  NAMES 

GRATIS 

00003 

HAFSVS 

00Q21 

INRT 

0O133 

FPGRAT 

0(1002 

EXTERNAL  SYHBqlS 

QllClOlOO 

01C03100 

oucosgT, 

PL'«T 

hercyr 

CHFelB 

00163  5YH80US 


lOENT 


ORATJE 
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SUBRSgTlNE  O^SVMaA8BLiXDlM|YDlM,NCHAR,JSU0> 

CATA(HH  ■ 0|C7»i  (h  I 0«08)i  (H  n 0,14)i  (8  « 0,0<)i 

*NC  » 3,0 

IQ0  » XASSr {INCHAR >-jSUB 

G8  TO{10,20i'!0i10|2Q)|O8 
ig  ANC  P ANC»ifC  « « «» 

20  CAI.U  5YMBeL(XDIMiYDIM,H,LA0BL»fl|O»8) 

X2  » XDIH*W»<NC#S*{ANC"Ii0>*0ji03 
YH  t YDIMyH 
IF{NCHAR,0T,C)Q9  T9  30 
X3  « X2 
00  TO  <0 

30  CALU  $YHBet.(X2iVH|HHilH0i0.0iU 
!F(NCHAR,1,T,5)Q0  T0  5q 
X3  ^ X2t3,0»{H*S) 

*-Q  CALL  $YHB0L{A3iYHiHHilH/i0i0iU 
5>0  RETURN 
END 


31300 

3140C 

31500 

31600 

31700 

31800 

31900 


32200 

32300 

32400 

32500 

32600 

32700 

328go 

32900 

33000 

33100 

33200 


lOENT  CHSyH 

PROCRAH  LENOTH  00174 

entry  P0INTS  0H5YH  OOOU 
EXYERNAt  SYMB018 

oaCDtcti 

SYNBfll, 

OC063  SYHSeiS 
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®’^PS®Ii<ft.PTO,PLPTH,BCO.N0CH,,vnj 

J ^ ^■-'59),  (130U.62B), 

C1MEN510N  LTMIq) 

NEWS  n XABSr^NH) 

LPTD  g ABSP(ftPTC) 

I.PTM  g FLPTM 
IF(NE,Gt,0)Q6  to  ?Q 
!F(l,H7M,NE.O)riO  TO  50 
NOCH  g 5 
ICAR  I lpTD 
GO  TO  «0 
90  NOCH  a 8 

ICAR  g UpTD*JOOO*l,RTH 
00  lF<l.PTO,LT,iCg,NeCH  * NOCH-l 
, D0  5 UZsIiNOCh 
5 I.TR(t.2)  > 0 
ACHAR  9 N8CH»2 

CALI.  0ULA8(rCAR|NCHAfljLTR<l)) 

KK  ? NOCH.l 
DO  10  JjaXK,Z 
10  LTR(JJ)  * IBLNK 

|F(NOCH,GT,6)LTR(N0CH*4)ajBLNK 
CO  TO(20,30>AEWS 
20  '<0«P  a lEAS 

l''<ftPTD,LT,C,OjK0Mp  I JHE5 
00  TO  4q 
30  KOHp  g ^joR 

JF <FLPTd,LT,C ,O)K0MP  I iSOU 
40  LTR{N0CH)  . KOMP  " ■ * 

CJ^jj^|jP*CK{NeCKATB|l),8CCl 

END 


PROGRAM  length 
EA'fRy  P0{f/rs  aoccoN 
CXTERML  SYM90L8  ^ 
OKIOIOO 
OOCOICT, 
RUlAB 

00J20 


00276 

0CQ22 


*0ENT  BCOCfN 


SUBROyTlNE  CENTR(MCiOL0C,OTLL) 

CgMMefj/CHfLOL/Ml.REs,  Hn^ER(50) ,HAXmEr,  HI,pAR, 
} naJ(Par,  n5,  MNC(5o»,  NTd, 

< NY(l-O),  NRY,  XL(50)i  YK(50>, 

3 NTAXISi  *'• 

LTA 


MmPAR<90), 

nLYi 

HXAXlSi 

YMB, 


CATAtHiO  S 0,08h  {SP  g 0,04» 

CTLL  • {mIOASP)RMC»8P 
CLOC^rnLL*,? 

RETURN 

END 


15000 

15100 

15200 

15300 

15400 

l55oo 

15600 
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PR0GRA>'  LC-NQTH 
Entry  points 
block  names 

EXTERNAL  symbols' 


00092  symbols 


CeNTR 

CHTL3L 

S 

0^004X00 

oacotcT, 


OOOS4 

00009 

Q0473 


IDENT 


CEnTR 


subroutine  oi-lasel 

C0Mm3N''m*»'3Y*/  C0RR<4)i 

1 iPRi 

2 YMAXi 

a RAH] I 

common/inpt/alatno,  alatnh, 

1 Al.GGNO;  AL0GNH, 

2 IN0TUSE<79) 
C0MM0IN/CHTL8L/H1NMER', 


RR« 

XMIN, 

HP, 

RAWI 

ALATNSi 

ALOGNSi 


DE, 

YMIN, 

ICRi 


ALATXO, 

AL0GXO, 


ALATXM, 

alogxm, 


ACHT, 

XHAXi 

lA. 

ALATXS, 

alogxs, 


C * . 1 1 


c 

C, 


CATA<N0  • 8)i 


MAXPAR'i 

NY(So)i 

XYAXISi 

LYN 

fKH  S 


MNMERI50)i 

nc, 

NRYi 

lYNi 


MAXHERi 

NNC19n>, 

XLISqJi 

Nl, 


MINPAR, 

NTC, 

YKISqJ, 

XTIB, 


MNPAR{9oJ, 

NLYi 

KXAXiSi 

YTiB’i 


0i07) 


C 

C 1 1 1 1 
100 


'•iio 


1^0 

iid 


ENEOEg  g 3,1419926535/190,0 
MIlN  3 »1 

|P(OELiLT,qnEDEG)nMIn  • 1 
LNG  « NH1N*1 
LTD  • NmIN*2 

call  OCOC6NOAL0ONO,AL0QNH,HlNMER,NC,LNGi 
CALL  8CDC6NtAL0GXOULeGxM»MAXMER,NTC#LNQ) 
CALL  0CDC'iNlALATNO,ALATNM,MJ,VPAR,NLYiLTOJ 
call  eCDC9fUALATX0,ALATxM,MAXPAR,NRYiLTD) 

set  up  minimum  merIoUn  label 

Yl  • YM1n»0,'209 
KNT  » O 

CALL  CENTR(NC,HTLl.TLUJ 
XI  g XM|N»HTLL 
K{  !l  O 
KNT  • KNTAI 

CALL  CHSYM(MlNMER,Xi,Yl,NC,3) 

SET  UP  INTBRHEPIATE  MErIOIAN  LABELS 
K|  » KUl 

KN  * XABSriNNCIKl )) 

call  centr<hn,htll,tll) 

Xl  » XL(KJ)4^TLL  , 

CALL  CHSYM{MNMeR(Kt)iXl|Yl,NNC<X|)»l) 
IF(Kl,LT,KXAXl3jCe  TB  120 

SET  UP  MAXIMUM  MERJOJAn  LABEL 
CALL  CENTR(NTC*HTLL,TLL) 

XI  « XMAXiHTLL 

CALL  CHSYM(MAXKER,Xl#yl,NTC|3} 

G9  TBaSoilflcJXNT 


26 


22600 

22700 

22800 

22900 

23ooO 

23100 

23200 

23300 
23400 
235oo 
236oo 
23700 
238oo 
23900 
24j 
24j 
24200 
24300 
24400 
24500 
24600 
2A7oo 
24800 
249qo 
25000 
29100 
25200 
29300 
25400 
25900 
25600 
29700 
29800 
29900 
26000 
26100 
26200 
26300 

i»Si 

26600 

26700 

26800 

26900 

27ooo 

27100 

27200 
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C 

C 


c 

c 

c 


ft.  SET  UF  FOR  IIPPER  B0f»OERi 
i60  Y1  » yHAX*0,C55 
G6  T9  100 


mm  set  up  reRIFlMHUM  PaRALLEI.  UA88L 
FIND  XI  PER  t^FT  6e«0ER 
leo  CALL  CEnTRCNUiHTlLiTLLI 
RNT  f 0 


200 


C 

C I . I , 
220 


290 


240 

c M M 
280 


290 

300 

C 

c » < M 

320 

390 


KY  1 0 
KNT  • KNT*1 

CALL  CHsyK(HJNPAR,XliYl|NLY|3) 

SET  UP  INTEFMBOIATE  PARALLEL  LABELS 

n Y • KY^l 
VI  A YK(KY)2HH 
ce  T8(230|24C)KNT 

find  XI  PER  left  eeROER 
KN  A xA0sr(NV(i<yn 
CALL  CENTR(K^ iHTLL.TLL) 

Xl  A XMIN«(TLL+0i125> 

CALL  CHSVM<HAPAR(Ky},Xl,Yl,NY(KY)il| 
IF(KY,LT,KYAXI3)68  T8  2Jo 

set  up  MAXJFUH  PARALLEL!  LABEL 
VI  A YMAX»fl,-2l 
00  T0{29O,3OC)KNT 

FIND  XI  FOR  left  80ROER 
CALL  CENTR(NPY,H7LL|TLL) 

5<1  3 XMlN^dLLAOiOSS) 

CALL  CHSyP.{HAXPAR,xl,yl|NRY,3) 

69  T9(32oi35cJKNT 

SET  UP  FOR  RIQHT  EOROER 
XI  A XMAX*0,C55 
09  T9  200 

return 

END 


27300 

27400 

27500 

27600 

27700 

27000 

27900 

28000 

2*100 


28200 

28300 

28400 

28500 

28600 


28700 

28800 

28900 

29090 

27100 


2’200 

29390 

29400 


295flo 

29600 


297oo 

29800 

29900 

30000 

30190 


30400 

30590 

3o6oo 

3o7oo 

30800 

30900 

31000 

31100 

31200 


PHSORAN  LENBTH 
ENTRY  POINTS  CHLA0EL 
eueCK  NAMES 

HAFSYS 

INFT 


„ , CHTLQL 

EXYERAAL  3YMB0LS 


oacoicT, 

B0IQ9N 

Cfi.'TR 


chjym 

OOHl  SVH86LS 


(10911 

00005 


lOENT 


0002X 

gOl33 

00473 


chlasel 
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BLODGETT  AND  MASSINGILL 


SUBR0UTINE  APACKINChS.STRHOiOUTLAB} 
CIMENSI0N  StFING<NnHS)»3UTLAaa) 
ENCQDE(a.iQliaUTLABl(STRlNQ(l)iI>tiNCHS) 
RETURN 

iOl  reRMATCBRi) 

BNO 


PR0GRAR  LENQ-^H 
ENTRY  P6INTS  AMCK 
EXTERNAL  SYliaOLS 

thenb, 

QdCOICT I 
ENC, 

aNSlNCL, 

00039  SYN8BLS 


pogoo 


I DENT  .APACK 


SU8R0UTINE  O>'F«L0(C0RA3»AU0C,NTRy) 
C3MM0N/|n/’T/AlATND,  0UHMY(5)i 


Al8QND» 

C0SM0N/CHHQt/H|AHE^  MN'mERjSO) 

1 MAXPARi  NCi  NNCe90)| 

2 NY<90>»  NRYi 

3 KYAXIS<  JYN# 

i LYN 

CIMENS10N  LTT<10) 

CATA(IBLNK  I 60600060606060900) 

KO0  R XASSFCN'TRY) 
oe  T0(io<^Oi^7O)KGe 
10  K3R0LO  >1  A0SF(A1.6GNO) 

jr<C0RAD,LT,e,O)KeR0LD  ■ KBR0liD«l 

IK  8 X 

KRT  8 NTRY 

NTRY  I 3 

lAX  * 1 

KXAXI5  8 0 

$ 9 7 9 3 0 ^ 

20  K0R0I.P  ? ABSMALATND)  , 

IF(C5RAD,LT|CiO)K0R0l.O  « K0RBLO*i 
IK  8 ; 

NTRY  5 3 
lAX  8 2 
KYAXl$  9 0 

30  call  RT0DM<O<RAPiK0RDOiC6RMN) 

K0RMN  ■ CSRMN 

K0RDO  8 XA95FIK0ROO)  , 

inK0RDQ,E(I|K6R0LOtARO'X®*^MN«N8|O)OB  ft  AO 
K K f ^ 

IF(K0RD5,LTrlOG)M<  I NK«1 
KOROID  8 KORPO 
IF(KRTilT,O)C0  T;  39 
ICAR  8 K0RnO»lflOO*K0RMN 
IF(C0RAD,lTrCiO)R9R0l’D  ■ K0ROO*l 
00  T0  50 
39  ICAR  8 K0RDQ 
NK  I NK«3 
00  T0  50 


DUhEECBA) 

MINPARi  MNPAR(90)i 
NTC»  NUY| 

YKl5n>i  KXAXiSi 
XYIP#  YTIB» 
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40 

90 

«0 


120 

190 

190 

160 

200 

290 


260 


300 


NK  3 


ICAR  a kORHN 

JK  f XA0SF(NK} 

Ce  60  UZaltJK 

LTT(LZ)  a I, 

CAU,  8ULA8|pCAR,IJX,LTTll)| 

00  T0|3ooil2Cil3Oi3oOil5oil6o)JK 

■ *8LNK 

G9  To  ZQO 

LTT(4)  f tTTC9)  » LTTJ6}  » IBLNK 
00  T0  200 

LTT(9J  ^ UTT(6)  • IBLnK 
00  T0  200 


LTT<4)  ^ laUA'K 

CALL  APACK(JH,LTT<1),ICaT> 

00  re(290|26C}lAX 


HNMERUK)  n ICAT 
NNC{IX)  P NK 
XLUK)  a AL0C 
KXAXIS  8 KXAXI$6l 
IK  a IK«1 
00  70  300 


KNPAR(IK)  a ICAT 

KY<IKj  a NK 
YK(|K)  a AL8C 
KVAXIS  a KVAXI961 
IK  a (K«1 
RETURN 
SNO 


19700 

19800 

19900 

20000 

2q3oo 

20A00 

20500 

20600 

20700 

20800 

20900 

21000 

21100 

21200 


21300 

21400 

21500 

21600 

21700 

21800 

21900 

22000 

22100 


22200 

22300 

22400 

22900 


PROGRAM  length 

Entry  points 

8L0CK  NAMES 


CHF0LB 


INFT 

external  5TMB0L5^^'^^^ 


OlCiOitQQ 

O0CDICT 

RTSOM 


9ULA9 

APACK 

00J45  symbols 


Q0340 

00016 


IDENT 


00133 

00473 


CHreLB 


subroutine  PT0RtC6O,AMlN,SeC,RAO) 

C5NVERSJ0N  67  OEGREES,  MINUTES,  SEC0NOS  1NT9  RADIANS 

PI*3, 141592654 

CEnD3SJGNF(f£G4AK IN/60, ♦SEC/3600, >#DEO> 

RAD?PI.DBOD/i90iO 

RETURN 

END 


59800 

59900 


56200 


563rtfl 

96400 


56500 
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PR6CRAH  LENQTH 
entry  P0INTS 
EXTERNAL  SYMBetS 


DTIR 

08CDICT. 


lOENT 

OOiU 

Q0003 


DT8R 


30 


subroutine  FlNnA{RR,X,Y,ALTX*.ALTN#Al.Gx.AUON*Ait^ 

this  SUORCUTINE  ID  USED  TO  CALCUL'mE  THE  SEHl-MAUSR  AXIS  9F  THfc 

shherrid  In  plotter  units  -n  no  further  conversions  are  necessa 
necessary  to  output  units  T9  THE  PLOTTER 

IF(XUO»10,3C 

CALL*  MERCTR U » 1 1 0 » RRj  ALTX«  ALOXi DUMMY • ^LAMX) 

CALL  MERCTRUiL.OiRRi ALTNiALQNiDUMMY»FLAHN) 

AaABSriY/IFLAMXiFLAMNJ) 

ITR  • 2 
GO  T0  ^0 

A?AaSr<X/tAU(:X*ALGNn 
RESTORE  BASE  MERIDIAN  AND  PARALLEL  IN  MERCTR 

AD  CALL  MERCTRll,A,fiR',ALTN|ALON,U0BOjO0D3) 

FIND  ACTUAL  HEIGHT  AND  WIDTH  OF  CHART 
CALL  MeR(;TR(C,*iRR,ALTX,ALQX;,XMAX#YMAX) 

CALL  M6RCTR<c»*i*’'’'**l-Tf^»*'-°N»XHlN»Y»lll) 

X » XM^X-XMJN 

Y « ymax.ympn 
RETURN 
ITR  • •! 

return 

END 


90 


56600 

56700 

56800 

56900 

57000 

57l00 

97200 

57300 

57400 


97700 

57800 

57900 


58000 

58l00 

58200 


98300 

98450 

58900 

58600 


98700 

58800 

58900 


PR80RAM  LENQtH 
ENTRY  POINTS 
EXTERNAL  SYMBOLS 


PPNDA 

08CDICT 

MERCTR 


00105  symbols 
subroutine  mercat 

COHrtON/HAPSYE/ 


1 

9 


00246 

00009 


I 


C9RR(Oi 

IPR| 

YMAXi 

RAH|i 


IDENT 


RRi 
XMIN, 
MR  I 
RAH! 


riNDA 


DELi 

YHIN( 

ICRi 


ACHt, 

XMAX, 

u, 


COMMON/INPT/ALATNOi  ALATNM,  ALATNSi  ALATXOi  Aj.ATjJH, 

1 ‘LBGNOi  AL0ONM|  ALOGNSi  AL5GXD»  ALOGXM,  ALOQXSj 

2 INPC,  KOMP,  WIDTH,  HITE,  DEBDVD,  DEGDVM, 

3 'teODVS.  NTIC,  TICLOTt  A1.TMOI5) , ALTMM{9  ALTMS  (5) , 

4 ALNMC;>(),ALNMM<8)iALNMS(5j,ALTS0(9),ALTSM|9),ALTSS{5)i 

) ALNS!2(5JiALNSM<5)iALNSS{5),ar{5)#  VR(5) 


CALL  riNDA<R8»HlCTHiHlTE«C8Rn<2JiC0RRUV»C0RR(4>iC3RR<3JiACHTilRf) 

iniRm30#i''Oii0o 


87800 

87900 

88000 

88100 

88200 

883oo 

88400 

88500 


88650 

887oo 

88800 

889oo 


89100 

89200 


30 


NRL  REPORT  7930 


100 


120 

130 

101 


CONTINUE 

RAHJ  K (C6Ri^(2>t»CeBR<l)J/HlTG 
BAWI  R {C0RR<4)sC8RR<3)>/HlOTH 

CALL  MERCTR(C,ACHT,Rft,CORRa),C0RR(3),XHlN,YMlN> 

CALL  MERCTR(C,AOl-T',RRiC0RR«2)iC0RR<4)iXMAX,yHAX> 

RETURN 

«RIT6(IPR,101) 

00  T0  12o 

PGRMATt*  JNSUFFICJENT  1NF0RHATI0N  T0  COMPLETE  PLAT*) 
END 


89400 

89900 

69600 

89700 

89S00 

89900 

90000 

90100 

90300 


PR8GRAH  LENGTH 

00071 

ENTRY  P0INTS 

MERCAT 

00014 

BL8CK  NAMES 

MAPSYS 

00021 

INPT 

00133 

EXTERNAL  SYMBOLS 

thend, 

oacotcT, 

FJiNDA 

merctr 

STN, 

lOENT 


MERCAT 


00104  symbols 

subroutine  MERCTR<1TiAiRR|AI.AT>AL0N|5{|Y^ 

1F<IT)10,20»10 
10  QM  n AL0N 
BP  R ALAT 
X8  9 A*9M 

E « S0RTP{1»C«RR*RR) 

E2  R E/2,0  . 

cINBP  ■ 



X n XB 

Y » Y0 
RETURN 

20  SINLA  « SiNFtALAT) 
eSlNLA  9 ErSINLA 
XA  R A'ALSN 

Y R YA»Y0 
RETURN 
END 
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90400 

90900 


9070C 
90800 
90900 
91000 
91100 
91200 
9l300 
9l400 
91900 
91600 
91700 
91800 
91900 
92000 
92100 
92200 
92300 
92400 


I 


BLODGETT  AND  MASSINGILL 


PReGRAH  length 
ENTRY  R01NTS  HBRCTR 
EXTERNAL  3YM00LS 

02C07m 

oacoiCT 

SORTr 

srNf 

L0cr 

cGsr 

00079  symbols 


Q02S1 

Q0009 


IDENT 


herctr 


SUBROUTlNe  TICX 


• T.MS  SUBROUTINE  ST8RES  uAlUES  IN  AN  ARRAY  THAT  ARE  USEO  POR  PLOTTING 

* TICK  MARKS  ON  .A  HERCATOR  0R  INVERSE  MERCATOR  PROJECTION  CHART 


CQMM0N/MARSYS/ 


cohmon/inpt/.alatnd 
aloc 

(•PC, 


C0RR(4)i 

RRi 

DELi 

ACHT, 

IPR| 

XMIN, 

YMIN, 

XMAX, 

YMAXi 

MR» 

ICRi 

14, 

RAHli 

RAHI 

1.  ALATNHi 

ALATNS* 

ALATXDi  ALATXM* 

ALATXS, 

I',  AL0GNM, 

ALOGNSi 

AL06XO,  AL0GXH, 

ALOGXS, 

K0MP| 

WIDTH, 

HITE,  DEGDVD, 

DEQDVH, 

i,  NTIC, 

TICLGT* 

ALTMD{9),ALTMM(3) 

,ALTMS<5 

ALNMC<5)iALNHM<5),ALNMSl5j,ALTSD(5j|ALTSM(5),ALTSS{5‘), 
AlNSC<9)Ui.NSM{5)»ALNSS(5j,aRI9),  VRI5) 


C0MMQN/riXORAT/DCL,NT|X 
C 

DEQH«CORR<4>tCQRR(9) 

XtNCR!XMAX«DCL/('CEQH*NT2X) 

C 

C«f*«TlCK  MARKS  ON  LOHER  OORDER  OP  CHART 
C 

lF(TlCLGT,LEiQiO)TlCLOTRt,P# 

106  TIC  • TICLOT 
Y0!|  YNIN 
M ■ 0 
XX  8 XM|N 
XX8XX«X(NCR 
IF<XX  .QE,  XMAX)90«19 
M«M«1 

ir{M»5)2’,l',4 
TICi»8,o«TIC 
00  TO  3 
ir(M-lO)3,0^0 
T I CjPij  j ntT  IC 
H a 0 
00  TO  3 
TICIPTIC 

CALL  PL9T(XXiYlr3) 

YA  ■ YOaTlCj 
CALL  PLeTCXXiYA,!) 

00  TO  lo 


00 

10 
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140000 

140100 

140200 

140300 

140400 

140900 

l4o6oo 

I4o7oo 

140«00 

140900 

141000 

141100 

I4l200 

l4i3oo 

141400 

141900 

1A1$00 

l4l7oo 

i-^iSoo 

141900 


14200 

14210 


142300 

142400 

142900 

142600 

142700 

142800 

142900 

143000 

143100 

143200 

l433oo 

143400 

143900 

143600 

143700 

143BQ0 

143900 

144000 

144100 

144300 

144900 


I 

1 


J 


n ntycy  r»or>r> 
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C**«*  TICK  marks  0N  upper  0eR[)ER  Qp  CHART 

V 

« YelyMAX®^'  f^MAX«0,Q005))55,35 

TICa.TlC’ 

GO  T0  60 




PRfiGRAM  UENQTH 
Entry  pqints 

eUGCK  NAHgS 


EXTERNAI,  SYMB0US 


HAPSYS 

INPT 

PrXGRAT 

Ql’COA^OO 

OaCDlCT, 

puer 


Q0135 

00003 

00021 

08133 

goooa 


I DENT 


T 

0012A  symbols 

SUBROUTINE  TlCY 


! tiiJk  '^SCD  PQR  PLOTTING 

• TJCK  H RKS  0N  4 WGRCATBR  0R  INVERSE  MERCATOR  PROJECTION  i;haRT| 


common/mapsys/ 

■ 


C0RR<A)i 

IPR, 

YMAx# 

RAHI, 


RR# 
XMIN, 
HR  I 
RAW! 


DEL. 

YMIN, 

ICR, 


ACHT, 

XHAX, 

lA, 


C0MM0N/INPT/ALA7NO;,  alatnm,  alatns,  ALATXD,  ALATXM,  alatxs, 

a AU0GXD,  ALOGXM,  ALOGXB, 

1 - R"MP,  WIDTH,  HITE,  DEnoVD,  DEGOVN, 

j n®*  NTIC,  TICUGT,  AUTMD (5 ) , ALTMMj? ) , al TMS (5 J 1 

a (5 ! iALNHM  (5  ) ,AlNMS(5)  , ALTSD(5)i,ALfSM(3),AL7^Sl5)i 

9 AlNSCJ5J.ALNSM(5),ALNSS(5i';ar<5>,  VR(5»‘  » IR  » , 


ALATXM, 

ALOGXM, 

DEGDVD, 


ALATXS, 

AL0GX8, 

DEGnyw, 


6C  ALAT?g0RR{l) 

YiNCRpDEL/NTiC 

ir{TlCLOT,LE.0,0)TICLQT|0,06 
166  TIC  • TICLOT 
XX  • XMIN 
Ml  ■ 3 

200  N « 0 

YY  » TMIN 

too  ALAT»ALATAYpNCR 


168000 
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20 

N=N*1 

l489go 

ir(N-5)5,4,l5 

149000 

4 

TICi»2jo.TIC 
GO  T9  « 

l49ioo 

149200 

15 

lF{N«lO)5il7tX7 

149300 

17 

TlCl>3,o*7IC 

l494flo 

N » 0 

149500 

GO  TO  4 

149600 

9 

T1C14TIC 

149700 

6 

CALL  PL0T<XX|YY,3) 

XA  a XX*T1CI 
CAUL  PL0T{XA,YY|2) 

149900 

35 

GO  TO  loo 

l50100 

99 

I»(XX  ,GE,  (XMAX»0,0009m03f48 

150200 

48 

XX»XMAX 

I5o3oo 

TIC»«TIC 

150400 

ALAT«C0RR{1) 

150500 

Ml  f 4 

150600 

GO  TO  200 

150700 

105 

CALL  PL0T{XMIN»YM;n,3J 

4q1 

RETURN 

154300 

END 

154400 

PR3CRAK  tENQTH 
entry  PQINTS  TPCY 
B1.3CK  NAMES 

MAPSYS 
INPT 

EXTERNAL  SYMBOLS 

Ol'COStOO 
Q8C01CT, 

MERCTR 
PL'5T 

0QiZ7  SYMBOLS 

subroutine  BLIABUNLABiKJiITEMP) 
riMENalON  ITEHPJKJ) 

ILA0.1NLA3 
C9  10  iKcl.M* 

K«KJ»IK*1 

ITE,MP(K)=M0DULA0,iO) 

il*bmlab/io 

IF<lLAO)i0,li|10 
10  CONTINUE 
u return 

END 


B0147 

QOOOS 

00021 

00133 


IDENT 


TICY 


13300 

13400 

13500 

13600 

13700 

13800 

I'^’OO 

lAOOO 

1<U00 

l'!200 

14300 


34 


NKL  HKPOKT  7<);!0 


lenqth 

t:\THV  PSiNTg 
EXTfcRiNH  Sy^^80l,S 


I'yiAB 


o'JCotgT, 

00033  syM00i.§^''^^^ 


00103 

00  003 


lOfcNT 


JJOD  I 180.0/Pt 
tSLTA  s O.OOCOOOl 

END 


BULA9 


128400 

128500 

128800 

128700 

128800 


120900 

129000 

129100 

129200 


rnovHAI’  (.ENOTH 

Entry  Roj^Tg  RTenw 
feXTERNit  symbols  ° 

Ol'ClOlOO 

OK03100 

00C33 


UU114 

OOQ03 


C 

C 


PJ=*3  1415926J3 
XRA[)iOfOl745029 
TEMW^aNQ 

(yr.xA)  1,10,1 
i ;j°'‘];^T.YAr/<XpT.XA> 

IE  (YRj.ya)  3,4,9 

^^NG?UtAN(Afl«,,/XRAo 

IE  <XPT,XA)  !ilO»iO 
^■'<GtE17XRAD 
CO  TO  14 

]E  (YPT.yA)  •7,4,8 

10  jr  (YPTsYA)  11,13,12 

11  *NG.27o,o  ***'»»iZ 

GO  TO  U 

12  AMO,9o,p 

60  TO  14 

13  ANG»»5M|» 

1<  RF'iUON 

BHD 


Q 3 
0 A 
0 7 

0 S 
0 9 

0 10 
0 11 
G 12 
0 13 
0 14 
0 IS 
0 18 
0 17 


0 20 
Q 21 
Q 22 
0 23 

Q 24 

a 23 

Q 28 
Q 27 

a 20 

0 29 

Q 30 
0 31- 
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IDENT 

PRGGRAH  LENQTH  00242 

ENTRY  points  ANCANO  00Q03 
EXIERNAU  SYMB01.J 

O0CDICT, 

ATANf 

00976  symbols 


ANQANO 


K 


! 

i 


I 

i 

I 

I 


I 


subroutine  riXANe(XPT,YrT,ANQ#CHAR,K9UNT.lLlNE,IPT’,  IUP) 

C annotates  P|X  numbers  I 2 

C t 4 

AHT,0,07 
NnG»2 

XRAD«0i01745329  I 9 

CNGnANG.90,  1 6 

If  (H0D(|<CUNT,2),NE:OJ  Q0  TO  1 17 

C ANNOTATES  PpGURBS  TO  T^E.  RIGHT  HAND  SIDE  0F  THE  DATA  POINT  I « 

XP»XPT  I 9 

YP»YPT  1 10 

GO  TB  4 I H 

C ANNOTATES  PfOURES  TO  THEi  UEPT  HAND  SIDE  Of  THE  DATA  POINT  I 12 

1 eN65{CNG*180i)«XPAO  I 19 

IP  (NPIG.CTiM  08  to  2 I lA 

HYPsAhT*(6,/7 , )«(NPIG#1)  I 19 

GO  TO  3 I i« 

2 HYPaAHT«{6,/7,)*NPlQ  I 17 

3 YY«HYP*S!NP(ENQ|  I iB 

XX*HYp»C0SPfGNQ|  I 19 

YPiyPT*YY  I 20 

XP.XPT*XX  I 21 


C ILINEsO  every  point  IUNE»1  EVERY  N TH  POINT  ANNOTATED 

c iUNE»2  Plot  value  series  no  annotation 

4 IPdLlNE. Die, 11,9 

11  IPIMODiKOUNTi  IPT),NE,1)  GO  TO  9 
CALL  PL0TIXP1,YP7,1UP) 

CALL  NUHBEflpP,YP,AHT,CHAR,CNG,2HP9l 
CALL  PL0TIXPT,YPT,3) 

CALL  symbol (APT, YPT,o,lo5iO»OtOi«l> 

GO  TO  7 

9 CALL  PL0T|XPT,YPT,1UPI 
GO  TO  7 

10  CALL  NUMBER  {XP,YP',AHTiCHAR,CNG,2HF9} 

CALL  PL0T(XPT,YPT,3) 

2b  CALL  SYMB8L{XPT,YPT,0,020»3i0i0»»l) 

GO  TO  7 


7 

CALL  PLOT  <XFT.YPT,3) 

I 

51 

IUP»2 

6 

RETURN 

1 

43 

END 

I 

44 

I 


i 

j 
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•—  V 


NRL  REPORT  7930 


LOiD 


PR0CR4P  LENQTh 
EN'TRy  P01NT5  PfXiWfl 

external  symbols 

ODC04:i0( 

(58C0ICT, 

PU5T 

^U^BER 

3Y^80L 

XMSO^ 

spNr 

C0SP 

0Q133  symbols 


00332 

00005 


IDEMT 


FIXAN0 


AB6RT,  'Uly  ON[T 


LUN  =«  «9 


ON  U'JASSlCNgn 
DfcCLARATjOMi  YOU  WAVE  Tfi?£ 

HAS  MOT  oeew  writtem  on.  i 

FOR  X @PTl8N, 


LUNp69, check  FORTRAN  CARO 
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BLODGErr  AND  MASSINGILL 


O »- 
O M 

n ^ 


CM  • 
O X 
o or 

Ui 


XXX 

XXX 

X a X 
XXX 
XXX 
XXX 
9 9 X 
XXX 
XXX 
XXX 
1X9 
XXX 
XXX 
XXX 
XXX 
XXX 
XXX 
XXX 
9 T X 
XXX 
XXX 
XXX 
XXX 
XXX 
XXX 
XXX 
V)  9 X X 
ox  X X 
X X X X 
«DX  X X 
<J  1 X X 
UIX  X X 
MX  X X 
^X  X X 
^9  3 X 
^X  X X 
— X X X 
X X X X 
XXX 
« X X X 

r«<  X X X 

«4X  X X 
»X  X X 
in  X X X 
XXX 
XXX 
XXX 
a X X X 
ox  X X 
XXX 
•XXX 
UX  X X 
UJ  X X X 
MX  X X 
9 X X 
»X  X X 

«n  a X X 

XXX 
-9  T X 
ZX  X X 
•-X  X X 
X X X X 
9 X X 
ox  X X 
ox  X X 
X z r 

•XXX 
MX  X X 
C Z X X 
ZX  X X 
XXX 
OZ  X > 
ox  X X 
Z Z X 
• 9 9 X 

UIZ  X z 
X X X X 
••X  z z 
»•  X Z X 
XXX 
ox  z z 
UIX  Z X 
MX  X X 
X Z X X 
< X X X 
-•X  X X 
UiX  Z X 
XXX 
WZ  X X 
« X X X 
M X X 7 
« X X z 
wiX  T X 
r«x  X X 
XXX 
XXX 
XXX 
r X X 
•-  X X X 
<Z  Z X 
XXX 
09X9 
UJ  X X 9 
X X X 
« X * X 
XX  X X 
«-»X  X X 
X X X X 
U A.  « « 
UiX  X X 
X X X 
XXX 
•4X  X X 
^ X X X 
t4X  X X 
OfX  X X 
OZ  X X 
OZ  X X 
XXX 
«K  X X X 
UI  X X X 
UIX  X A 
X X X X 
^X  X X 
^ X » 3 
XXX 
UiX  X X 
O X X X 
<r  X X X 
u*  *.  X X 
:9X  X X 
OM  M M 
UIM  X X 
V>  n M 
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